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Abstract. We investigated the physical and thermal properties of 65 Cybele, one of the largest main-belt asteroids. 
Based on published and recently obtained thermal infrared observations, including ISO measurements, we derived 
through thermophysical modelling (FPM) a size of 302 x290 x232 km (± 4%) and an geometric visible albedo of 
0.050±0.005. Our model of a regolith covered surface with low thermal inertia and "default" roughness describes 
the wavelengths and phase angle dependent thermal aspects very well. Before/after opposition effect and beaming 
behaviour can be explained in that way. We found a constant emissivity of 0.9 at wavelengths up to about 100 /im 
and lower values towards the submillimetre range, indicating a grain size distribution dominated by 200 fim 
particle sizes. The spectroscopic analysis revealed an emissivity increase between 8.0 and 9.5 /jm. We compared 
this emissivity behaviour with the Christiansen features of carbonaceous chondrite meteorites, but a conclusive 
identification was not possible. A comparison between the Standard Thermal Model (STM) and the applied TPM 
clearly demonstrates the limitations and problems of the STM for the analysis of multi-epoch and -wavelengths 
observations. While the TPM produced a unique diameter/albedo solution, the calculated STM values varied by 
±30 % and showed clear trends with wavelength and phase angle. 65 Cybele can be considered as a nice textbook 
case for the thermophysical analysis of combined optical and thermal infrared observations. 
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1. Introduction 

The asteroid 65 Cybele is a main-belt asteroid and 
' the main representative of a group of asteroids with 
semi-major axis in the range 3.27<a<3.70 AU, eccen- 
tricities e<0.30 and inclinations of i<25° (Gradie et al. 
1989|. 65 Cybele is generally considered as one of the 
10 largest asteroids, with a published diameter between 
230 (Dobrovol'skij HHHOI Tavlor 11551)1 and about 330 km 
(Green et al. 1985). Although discovered in 1863, it took 
until 1980 before its rotational properties were determined 
as one of the last remaining large minor planets (Schober 
et al. I1980|l . Recent analysis of many lightcurve observa- 
tions lead to a rotation period of about 6 hours, a unique 
spin vector solution pair with a retrograde sense of rota- 
tion and an elongated body shape (EriksonpOOO). 

Send offprint requests to: T. G. Miiller 

* Partly based on observations with ISO, an ESA project 
with instruments funded by ESA Member States (especially 
the PI countries: France, Germany, the Netherlands and the 
United Kingdom) and with the participation of ISAS and 
NASA. 



Based on its low albedo and the flat to slightly red- 
dish spectrum in the 0.3 to 1.1 jum range, it was classi- 
fied as P-type object (Tholen 19c25|) and due to a weak 
ultraviolet absorption feature short ward of 0.4 iim as C0- 
type (Barucci et al. 11987)1 . Recently Bus l|1999[> classified 
65 Cybele as Xc type, based on reflectance spectra cover- 
ing the wavelength interval from 0.44 to 0.92 /im. These 
taxonomic types indicate hydrated silicates, carbon and 
organics on the asteroids surface (Gaffey et al.|l989). 

Another interesting aspect comes from the occulta- 
tion measurement: Dobrovol'skij l|1980)l reported a sec- 
ond event during the occultation of AGK3 +19 599 which 
was attributed to a satellite of 11km diameter, located 
917 km from the center of 65 Cybele (in the case of a cen- 
tral occultation). So far, no confirmation of the satellite 
through other techniques was published. It was one of the 
targets in the HST search for binaries 'Imaging Snapshots 
of Asteroids' Storrs et al. (1993) & Storrs (priv. comm.), 
but with unknown outcome so far. No radar measurements 
are available yet. 

Tedesco & Veeder $TW^ and Tedesco et al. plH^ 
derived from IRAS observations diameter and albedo val- 
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ues of 237.26±4.2km and 0.0706±0.003, respectively. The 
only direct size information from the occultation event in 
1979 resulted in a diameter of 230±16km, with a largest 
measured chord of 245 km (Dobrovol'skij 1980; Taylor 
119811 Xiu-Yi[Tn7H|). Both values, the radiometric and the 
occultation diameter, agree nicely. But ground-based N- 
band observations (Morrison & Chapman 119761 Morrison 
115771 Bowell et al. IT575I Green et al. QUI lead to di- 
ameters larger than 300 km. The published radiometric 
diameters and albedos were all based on the concept of 
the Standard Thermal Model (STM; Lebofsky et al. [TUM 
Lebofsky & Spencer 1989), but using different values for 
the infrared beaming parameter 77, the phase angle cor- 
rection and the emissivity. Using for example the STM 
77- value of 0.756 (Lebofsky et al. 1986) would lower some 
of the earlier 300 km diameter values to about 250 km, 
but the Green et al. (1985) values would still be close to 
300 km. The higher phase angle correction of Bowell et al. 
(1979) would again increase the diamter of all measure- 
ments taken at large phase angles. 

The diameter is the most fundamental parameter, but 
in case of 65 Cybele it has very large uncertainties, much 
larger than for other asteroids of comparable size. Where 
is this uncertainty coming from? Are the data sets too 
different or some of them unreliable? What role play the 
STM concept, the simplification in the shape, the missing 
thermal and rotation effects or the various other model 
input parameters? 

The first thermal spectrum of 65 Cybele was published 
by Green et al. ([1985 ) . They did not see any pronounced 
emission features in the 8 to 13 y«m wavelength range. In 
the years 1996-1998 65 Cybele was extensively observed 
by ISO with the scientific goal of providing new insights 
into the chemical composition and the mineralogy for a 
sample of main-belt asteroids. Are there any spectral fea- 
tures visible? If yes, what are the meteoritic and miner- 
alogic counterparts? 

With the availability of a recently established thermo- 
physical model (TPM) by Lagerros lTi59l)l 115571 IT55£jl we 
are now in a position to combine the information on the 
asteroid, its rotational behaviour, different observing and 
illumination geometries with the thermal measurements 
taken at different wavelengths and different times. 

In the following, we show the results of the re-analysis 
and re-calibration of the existing thermal observations 
from ground- and satellite-based instruments, including 
photometric and spectroscopic measurements (Sect. [2J - In 
Sect.|21we present the basic parameters of our TPM anal- 
ysis and derive radiometric diameters and albedos for all 
photometrically reliable data sets (Sect. l3~2"f . In the spec- 
troscopic analysis (Sect. l3~3*|l we focused on the modelling 
of the shape of the thermal continuum emission which is 
strongly connected to thermal effects on a regolith covered 
surface. The TPM analysis section is followed by exten- 
sive discussions of radiometric and spectroscopic results 
(Sect.EJ. We used a large sample of infrared observations, 
spanning a wide range of phase angles and wavelengths, to 



demonstrate the possibilities and limitations of the ther- 
mophysical modelling of asteroids (Sect.[SJ). 

2. Observations in the Thermal Infrared 

2.1. Re-analysis of previous measurements 
TRIAD Observations 

The original TRIAD 1 fluxes were based on three 
65 Cybele observations from March 18, 1975, 08:24 UT, 
with an average N-band magnitude of +0.85±0.05mag 
(Morrison & Chapman 1976). We used the Rieke et al. 
( 1985) absolute photometric system for the standard stars 
to re-calibrate this flux. This led to a new N-band value 
of +0.92 mag, corresponding to a colour corrected flux at 
10.6/im of 14.55±0.87Jy (measurement and calibration 
errors have been added quadratically). 

IR Spectroscopy from UCL 

Green et al. <| 198511 presented the first 8- to 13 /jm spec- 
trum of 65 Cybele. It was taken on April 18, 1982 and cal- 
ibrated against a Boo. We re-calibrated the 10.0 /im flux 
value to 25.5±2.9Jy with the Cohen et al. I|1995|l model 
for a Boo (the Green et al. value for a Boo is 15.0% 
higher than the Cohen et al. value at 10.0 /im). The re- 
calibration resulted in flux values of 14.4±2.0 Jy at 8.5 /urn, 
25.5±2.9 Jy at 10.0 /_tm and 50.9±4.4Jy at 12.0 /im for the 
given epoch. 

IRAS Observations 

IRAS (Beichman et al. 1988J) scanned 6 times over 
65 Cybele, each time with all 4 wavelength bands. The 
latest flux values and the corresponding analysis were pre- 
sented by Tedesco et al. H2002|l . We used the colour cor- 
rected monochromatic IRAS fluxes (without applying the 
STM-related band-to-band corrections) for our reanalysis 
(see also Miiller & Lagerros 1998). 

IRTF Observations 

We included in our analysis 3 IRTF measurements (D. 
Osip, priv. communications, see Table fj]): 

Table 1. IRTF observations of 65 Cybele from Dec. 19, 
1996 



Julian Date 


A c [/im] 


Flux [Jy] 


2450436.71813 


10.5 


7.0 


2450436.72042 


10.5 


6.7 


2450436.72350 


18.0 


17.8 



No observational or calibration errors were given, but 
the data were considered as a high accuracy dataset (Osip, 

1 Tucson revised Index of Asteroid Data 
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Fig. 1. PHT-S spectrum of 65 Cybele. The TPM predic- 
tion is plotted as dashed line (with -D e // = 296.5 km, 
Pv — 0.043 and default thermal parameters as described 
in Sect. l3.l|) . The error bars are the signal processing errors 
(between 2 and 6 %) added quadratically to the assumed 
calibration error of 5 % per pixel. The solid line in the 
lower part was produced through a 10 element smoothing 
of the PHT-S to TPM ratio (see Sect. El). 

priv. comm.). We assumed a 10 % flux uncertainty for the 
N-band and 15 % for the Q-band photometry. 

UKIRT Observations 

The UKIRT mid-IR N- and Q-band data of 65 Cybele 
have been obtained in August 1996 with the MPIA 
5-25 /mi camera MAX (Robberto & Herbst HM^ . 
Observations, data reduction and calibration are described 
in Miiller & Lagerros (192EJ). 

ISOPHOT photometric data 

65 Cybele was measured at wavelengths between 65 and 
200 /mi with ISOPHOT (Lemke et al. llTM)) at 3 different 
observing epochs. Observations, data reduction and cali- 
bration are discussed in Miiller & Lagerros (1998 2002 ). 

2.2. ISOPHT PHT-S Measurement 

The observational details of the PHT-S measurements are 
given in TableEl The PHT-S short wavelength part (2.47- 
4.87/mi) covers the transition between reflected light and 
thermal emission. It has been excluded in our thermophys- 
ical analysis. The data reduction of the 5.84-11.62 (im part 
was done in a standard way, using PIA 2 (for details of the 
PHT-S data reduction steps see Dotto et al. 2002J). There 
was no dedicated background measurement, we therefore 
used COBE-DIRBE measurement results (Hauser et al. 
1998| as replacement. The color corrected, weekly maps 
from the COBE-DIRBE data sets were adjusted for the 
PHT-S beam size. The two DIRBE values at 4.9 and 12 /xm 
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Fig. 2. CAM-CVF and CAM filter observations of 
65 Cybele. The TPM prediction is plotted as dashed line 
(with D ef f = 278.8 km, p v = 0.049 and default ther- 
mal parameters as described in Sect. 13. if) . The lower part 
shows the resulting CAM filter flux densities divided by 
the corresponding TPM prediction (see Sect. [3j) ■ 



were then connected using a 266 K blackbody (Abraham 
et al. I1999fl and subtracted from the measured spectrum 
of 65 Cybele. The background contribution was between 
4 and 10 % in the PHT-S wavelength range. 

Figure^shows the reduced and background subtracted 
result of the PHT-S observation together with the corre- 
sponding TPM prediction. Typical signal processing un- 
certainties were between 2 and 6%. The absolute calibra- 
tion error for bright point-like sources in staring mode is 
less than 10% (Laureijs et al. 2002). For our analysis we 
added quadratically the individual signal processing errors 
and an absolute calibration error of 5%. The lower part 
of Fig. n shows the ratio between the observed spectrum 
and the TPM prediction (see Sect.^J). 

2.3. ISOCAM CVF & Filter Observations 

On Christmas day 1997, a for- and backward scan was 
performed in the 5 to 17 /im wavelength range with the 
ISOCAM Circular Variable Filter (CVF) (Cesarsky et al. 

1996) . The scan was performed skipping one CVF position 
per CVF step, corresponding to steps of approximately 
0.2 /im for the 9-17/mi and 0.1 /im for the shorter wave- 
length range. The pixel field of view was 3" leading to a 
total field of view of 96". The data reduction was done 
within the CAM Interactive Analysis package (Ott et al. 

1997) and followed the same strategy as was done for the 
standard stars observations, used for the calibration of 
the Spectral Response Function (Blommaert et al. [2001 ) . 
The flux calibration was done using the OLP vlO SRF 3 . 
Aperture photometry was applied avoiding the contribu- 
tion from reflected light on the detector (Blommaert et al. 



ISOPHOT Interactive Analysis 



Off Line Processing version 10, Spectral Response Function 
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Table 2. PHT-S, CAM CVF and CAM filter observations of 65 Cybele. OTT is the On Target Time. 







Start 


A-range 


OTT 






Instrument 


TDT 


Observation 


H 


[sec] 


Field of View 


PHT-S 


37300324 


23-Nov-96 12:40:26 


5.84-11.62 


364 


24" 


x 24" 


CAM-CVF forward (down) 


77002401 


25-Dec-97 02:08:41 


17.00-4.96 


2996 


96" 


x 96" 


CAM-Filter 


77002402 


25-Dec-97 02:59:19 


6.0-14.9 


3118 


48" 


x 48" 


CAM-CVF backward (up) 


77002403 


25-Dec-97 03:52:01 


5.02-17.00 


2970 


96" 


x 96" 



2001 ) . The strategy of combining a for- and backward scan 
allows to check observed spectral features and also allows 
to a certain extent the correction of the ISOCAM transient 
response. We applied the Fouks-Schubert method (Coulais 
& Abergel E96j to correct the transient of the signal. 
The CVF scan starting at 5 /im after the fixed filter mea- 
surements suffers strongly from a remnant which is not 
corrected for by the transient correction method. For this 
reason we judge the data for 5-7 /im in the up-scan as 
less reliable. 

In the time between the two CVF scans nine fixed filter 
measurements were performed in the 6-15 (im wavelength 
range. The zodiacal background becomes important at 
wavelengths longer than 9 fim. Therefore, dedicated back- 
ground measurements were done in the long-wavelength 
filters. They were performed before 65 Cybele itself was 
observed to avoid remnants of the 65 Cybele measurement 
because of the transients. The on-source measurements 
were long enough to allow stabilisation of the signal so 
that no correction had to be applied (Blommaert et al. 
2000}. The obtained flux densities have been colour cor- 
rected using the average CVF as input spectrum. 

The resulting monochromatic flux densities are listed 
in Table I3J These values are also plotted together with 
error bars in Fig. [21 The lower part of the figure shows 
the ratio between the observed CAM filter values and the 
corresponding TPM prediction (see Sect.EJ). 

2.4. ESO-TIMMI2 N- and Q-band Observations 

We obtained additional N- and Q-band observations of 
65 Cybele using the TIMMI2 camera (Reimann et al. 
2000), mounted on the ESO-La Silla 3.6m telescope. 
Table 0] summarises the observational results 4 . 

We did a standard data reduction for the chopping- 
nodding imaging mode. Asteroid and calibration star mea- 
surements have been taken within 30 min, all Nl inte- 
gration time were 144 sec, the Ql integration times were 
108 sec. The relatively large uncertainties in table 0] are 
due to the variable atmospheric conditions. The flux er- 
rors are estimates based on the monitoring of the reference 
stars. No colour or airmass correction was performed, both 
effects are only on a 1-2 % level and much smaller than the 
photometric errors. 



4 Based on observations collected at the European Southern 
Observatory, Chile; ESO, No. 69.C-0090 



Table 3. ISOCAM filter measurements. The monochro- 
matic flux values are background subtracted and colour 
corrected. The 1-cr errors of the photometry are between 
1 and 4%. The colour correction terms (Col. CC) have 
been determined from the averaged CVF-spectrum. CVF 
uncertainties produce errors in the colour correction val- 
ues of about 5 % in the LW4 and LW2 filters and about 
1 % for the others. The uncertainty in flux conversion is 
approximately 5 %. All errors have been added quadrati- 
cally (Col. 'Error'). The measurements were taken on 25 
December 1997 between 3:00 and 3:50 UT. 



filter 


Xref [A*m] 


CC 


Flux [Jy] 


Error [Jy] 


LW4 


6.0 


1.08 


0.368 


0.028 


LW2 


6.7 


1.32 


0.964 


0.073 


LW5 


6.8 


1.00 


0.899 


0.052 


LW6 


7.7 


1.06 


2.27 


0.14 


LW7 


9.6 


1.03 


6.85 


0.40 


LW8 


11.3 


1.01 


12.74 


0.73 


LW10 


12.0 


0.91 


15.20 


0.82 


LW3 


14.3 


0.99 


23.49 


1.25 




14.3 


0.99 


23.39 


1.22 


LW9 


14.9 


1.00 


27.09 


1.42 




14.9 


1.00 


26.92 


1.50 



Table 4. Summary of the ESO-TIMMI2 N- and Q-band 
observations. All measurements were taken between air- 
mass 1.02 and 1.24. The calibration stars were HD 178345 
and HD 169916, with model spectra from Cohen et al. 



date 

(UT) 


time 


filter 


Are/ 

[Mm] 


Flux 

[Jy] 


Error 

[Jy] 


10- Jun-02 

11- Jun-02 


07:42 
08:47 
08:51 
10:02 


Ql 
Ql 
Ql 
Nl 


17.75 
17.75 
17.75 
8.60 


74.9 
57.4 
57.7 
8.8 


±?« 

5.0 
5.0 
0.5 



3. Thermophysical Model Analysis 

We applied the thermophysical model (TPM) developed 
by Lagerros l|1996l 119971 11998|l and extensively used and 
refined by Miiller & Lagerros fiMfl I2UB2) | . The TPM in- 
cludes shape, rotational behaviour and thermal behaviour 
of a given asteroid at the specific observing and illumi- 
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nation geometry. It allows to derive radiometric sizes and 
albedos from thermal infrared observations, taking shape 
information, geometry and thermal aspects into account. 

3.1. Model input values 

The H 5 and G 6 values are taken from Lagerkvist & 
Magnusson l|199U|) . with H = 6.70 ± 0.15 as the weighted 
mean value of four apparitions and G = 0.09 ± 0.05. 
Erikson (2000) presented a good quality solution pair for 
the spin vector based on a retrograde sense of rotation 
and a synodic rotation period of 6.1 hours (see Table EJ- 
However, due to noisy and under-sampled input data, the 
spin vector pair is still considered only as class II, indicat- 
ing that higher quality data might improve the situation. 

The sense of rotation can in principle be determined 
from the anlysis of thermal observations before and after 
opposition (e.g. Miiller 2002 ). In cases where two spin vec- 
tor solutions with different senses of rotation are given, the 
thermophysical analysis allows to determine the correct 
one. But to distinguish between two spin vector solutions 
with the same sense of rotation is very difficult on basis of 
infrared data. The before-/after opposition effect is then 
not very pronounced. Our TPM analysis for 65 Cybele 
gave almost identical diameter and albedo values for both 
spin vector solutions. For the first solution (first line in 
TableEJ we obtained a slightly smaller standard deviation 
around the mean diameter and albedo values. Therefore 
we accepted solution 1 for all following calculations. The 
corresponding ellipsoidal shape model solution is closely 
connected to this spin vector solution. 

The determination of thermophysical quantities such 
as thermal inertia, beaming parameters or emissivity re- 
quire large samples of infrared measurements (Miiller & 
Lagerros ll998fl . but for regolith covered, large main-belt 
asteroids, default thermal values seem to work fine (Miiller 
et al. I1999[) . They found a thermal inertia value of T = 
15 J m -2 s~ 0,5 K _1 for 1 Ceres which they also used suc- 
cessfully for other main-belt asteroids (Miiller & Lagerros 
2002}. We adopted this T- value, as well as the derived 
beaming model parameters / = 0.6 (fraction of surface 
covered by craters) and p — 0.7 (r.m.s. of the surface 
slopes). These beaming values correspond roughly to the 
"default" STM beaming parameter (Lebofsky et al. H986[) 
if the heat conduction is neglected (Lagerros[T929 . For the 
emissivity we used a constant value of e = 0.9 at all wave- 
lengths between the near-IR and 100 /im. At longer wave- 
lengths we used a wavelength dependent emissivity with 
decreasing values from 0.9 at 100 /im to 0.8 at 400 /im. 

3.2. Radiometric analysis 

Based on TPM parameters from Sect. 13.11 we derived ra- 
diometric effective diameter and albedo values from all 



5 Absolute magnitude in the H-G magnitude system 

6 V-band slope parameter in the H-G magnitude system 



photometric and spectroscopic measurements in Sect. [21 
The results are given in Table 

The typical observational and calibrational errors of 
individual measurements are between 5 and about 20% 
in flux density. The translation into diameter leads to 3 
to 10 % uncertainties for individual radiometric diameters. 
A 10 % change in radiometric diameter for 65 Cybele re- 
quires an albedo modification of 0.010 (a larger diameter 
requires a lower albedo and vice versa). The D e ff values 
per instrument are within ±8 % of the weighted mean, 
only the UKIRT value is significantly off for unknown rea- 
sons. The ISOPHOT data were taken under the same as- 
pect angles as the UKIRT data, but do not show such 
high D e ff values. The highest quality photometric data 
are the CAM filter measurements. The derived diameter 
and albedo values are very close to the weighted mean 
values of all 10 sets of independent observations. 

The effective diameters are connected to the shape 
model via the equal volume sphere formula D e ff = 
2(a6c) 1/3 (see also Miiller & Lagerros ll998jl . Shape, object 
orientation (represented by spin vector and aspect angle 7 ), 
illumination, rotation and thermal effects are taken into 
account in the calculation of D e ff and py- The calcu- 
lated weighted mean effective diameter of 273.0 km corre- 
sponds to an absolute ellipsoidal size of 2a x 26 x 2c = 
301. 9x290. 3x232. 2km for the given spin vector solution 
(first line in Table |3J. 

3.3. Spectroscopic Analysis 

The original idea of the spectroscopic observations was to 
get new insides into the chemical composition and min- 
eralogy of this main-belt asteroid. Mainly the nature of 
silicate features can be studied in the wavelength range 
between 8 and 12 /im (Dotto et al. 2000). Green et al. 
(1985) found no prominent emission features in the N- 
band spectrum of 65 Cybele. 

Dotto et al. l{2"UU2"jl and Barucci et al. l{2"Ucl2"jl anal- 
ysed the spectroscopic features by dividing the observed 
spectrum by the corresponding TPM prediction under 
the assumption of a perfect emissivity at all wavelength 
(e = 1.0). Here, we follow this procedure. Figure |31 shows 
the PHT-S and the CAM-CVF measured spectrum di- 
vided by the corresponding TPM with e = 1.0. Both spec- 
tra agree in spectral shape within the given observational 
and instrumental uncertainties. 

The PHT-S spectrum (Fig.nj was derived from 64 in- 
dividual pixels which have a completely independent cali- 
bration for each pixel. The PHT-S data are therefore not 
ideal for an analysis of very low level spectroscopic fea- 
tures. The CAM-CVF scans are much better suited for 
such an analysis. The source is always seen by the same 
pixels and with the variable filters the source spectrum 
is scanned. Despite the disadvantages of the PHT-S spec- 
trum for the analysis of shallow features, a 10 element 

7 Angle between the rotation axis of 65 Cybele and the radius 
vector to Earth 
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Table 5. Spin vector parameters of 65 Cybele (Erikson|2000). 



Asteroid 


Spin Vector 


Sidereal Period 


Absolute rotational phase 


Ellipsoidal model 




Ao Ad 


Psid 


7o 


a/b b/c 


65 Cybele 


56±15 -25±15 


0.2572901±0.0000004 


-38° at T = 2446420.5 


1.04±0.05 1.25±0.10 




237±15 -16±15 


0.2572901±0.0000004 


-36° at T = 2446420.5 


1.04±0.05 1.24±0.10 



Table 6. D e f f and pv values, together with the standard deviation of all calculated values per instrument (for TRIAD 
only the flux uncertainty has been taken into account). The values have been derived through the TPM with default 
input parameters for shape and thermal behaviour (see lfl.ljl . 



Instr. 




A 


Asp. 


D.ff 


Pv 


(# of Obs) 




H 


t°] 


[km] 




Photometry 


TRIAD (1) 




N-band 


70 


256.4± 7.5 


0.058±0.003 


IRAS (24) 




12,25,60,100 


68-71 


262.4±14.0 


0.056±0.006 


UKIRT (4) 




N- & Q-band 


87 


(329.1±19.9) 


(0.035±0.004) 


IRTF (3) 




N- & Q-band 


94 


281.5±15.7 


0.048±0.006 


PHT-P/C (13) 




60.. .200 


85-94 


269.2±10.0 


0.053±0.004 


CAM Fil (11) 




6.0... 14.9 


144 


278.8± 7.0 


0.049±0.003 


TIMMI2 (4) 




8.6, 17.75 


66 


260.5±19.2 


0.057±0.008 


Spectroscopy 


IRSPEC (3) 




8.5,10.0,12.0 


26 


264.8± 9.6 


0.054±0.004 


PHT-S (12) 


6.O.. 


.11.5 (0.5 /im steps) 


89 


296.5± 8.3 


0.043±0.002 
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273.0±11.9 
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Fig. 3. Comparison between PHT-S (diamonds) and 
CAM-CVF (dashed line with error bars) spectra divided 
by the corresponding TPM predictions with e — 1.0. The 
solid line was produced by a 10-element smoothing of the 
PHT-S emissivity values. The CAM-CVF is the average of 
the up and down scans divided by the model prediction. 



Fig. 4. The observation over TPM ratios for the photo- 
metric measurements plotted against the corresponding 
wavelengths (see Fig.Hfor STM ratios). The IRAS values 
are given without error bars for clarity. All TPM calcula- 
tions were done with the resulting weighted mean values 
from Table |SJ 



smoothed PHT-S spectrum resembles the CVF spectrum 
very well. 
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Fig. 5. The observation over TPM ratios for the photo- 
metric measurements plotted against the corresponding 
phase angles. All TPM calculations were done with the 
resulting weighted mean values from Table El 



4. Discussion 

4.1. Physical parameters 

Our reanalysis of 10 independent data sets (see Sect. 0) 
give now very consistent radiometric result through 
the TPM. We determined the effective diameter to 
D e// =273.0±12.0 and the albedo to p y =0.0 50±0.005. 
The discrepancy to the Tedesco et al. H2002fl results is 
mainly related to the spherical shape simplification of the 
STM. Thermal effects, the beaming model, phase angle 
corrections and band-to-band corrections also contribute 
to the deviations. The occultation diameter was derived 
from a circle fitted to 3 measured chords. The uncertain- 
ties of the measurements and the limited information from 
a 3-chord occultation event also allow for a fit with a 
larger slightly elongated ellipse, corresponding to the pro- 
jected cross section of our derived radiometric size values. 
The fact that different data sets, taken at different wave- 
lengths, different observing and illumination geometries 
give now one consistent solution is certainly in favor of 
our larger diameter. And, all IRAS observations fit nicely 
into that picture, including the 100 /xm data which were 
not used for the Tedesco et al. (2002) calculations. Here, 
we would like to emphasize that we did not apply any 
band-to-band corrections or any other wavelengths-, phase 
angle or beaming corrections, like it was done for the STM 
IRAS applications. 

Despite the usage of best available ellipsoidal shape 
model and the nice agreement on average between the 
TPM and the observations, we believe that the unknowns 
of the true shape limit the accuracy of TPM predictions 
for individual epochs. Deviations between the true object 
cross section and the projected model cross section are 
easily possible, but so far no improved shape model (e.g. 
from radar or direct imaging methods) is available. 



The shape model is only useful together with the spin 
vector solution. Erikson l|2000|l gives 2 solutions as result 
of lightcurve inversion methods (see Table EJ), but both 
solutions have the same sense of rotation. The determina- 
tion of the correct solution on basis of radiometric analysis 
is therefore very difficult (see Sect. I3.ll for details). Here, 
we accepted solution 1 for all calculations. 



4.2. Thermal parameters 

We used a flat emissivity model with e = 0.9 over all wave- 
lengths out to 100 /mi, at longer wavelengths the model 
emissivity is slowly decreasing towards values of 0.8 at sub- 
mm wavelength. Previously published emissivity values of 
1.0 at around 20 /mi were based on ground-based Q-band 
observations and IRAS 25 /mi fluxes together with very 
uncertain beaming values (Miiller & Lagerros ll998[l . With 
our modification of taking a constant 0.9 value at mid- 
IR wavelengths we investigate the emissivity behaviour in 
this critical region where beaming effects play an impor- 
tant role (Dotto et al. I2000J) . In Fig. 0] there is no trend 
with wavelength visible out to about 150 /mi. The constant 
'observation over TPM' ratios in Figs. ^ an( i ID confirm 
our assumption of a constant e = 0.9 at mid-IR wave- 
lengths. The 5 ISOPHOT measurements at wavelength 
beyond 150 /mi (calibrated against a few stars and the 
planets Uranus and Neptune) show fluxes which are about 
10% lower than the TPM predictions (Fig.^J). Miiller & 
Lagerros (20Q2J| see nice agreement with this current emis- 
sivity model for 1 Ceres, 2 Pallas, 4 Vesta and 10 Hygiea 
out to 200 /tm, although for 2 Pallas and 4 Vesta there were 
indications for a similar trend at the very longest wave- 
lengths. The observed far-IR fluxes are to a good approx- 
imation directly proportional to the hemispherical model 
emissivity (Miiller & Lagerros 1998). Instead of a slow de- 
crease from 0.9 at 100 /tm to 0.8 at 400 /tm, the 65 Cybele 
observations now indicate a more abrupt decrease in emis- 
sivity to values around 0.7-0.8 in the wavelength range 150 
to 200 /tm. 

Redman et al. (1992| IT39"8"j) found effective emissivi- 
ties between 0.5 and 0.9 in the sub-mm range for 7 large 
asteroids. They attribute differences in the emissivity be- 
haviour to different optical depth of the warm surface ma- 
terial and its density. Scattering processes by grains within 
the regolith reduce the emissivity in a wavelength depen- 
dent fashion. For the Moon the far-IR emission is also 
suppressed by scattering processes (Simpson et al. 1981). 
Scattering processes are most effective if wavelength and 
grain sizes have similar dimensions. Following this inter- 
pretation, our far-IR low emissivities indicate that the 
dominant particle size in the 65 Cybele's regolith (and 
maybe also for 2 Pallas and 4 Vesta) would be around 
200 /im. More observations are needed to identify the pre- 
cise emissivity behaviour. Well-calibrated spectroscopic 
observations in the far-IR and sub-millimetre (e.g. with 
HERSCHEL) will then allow a detailed description of par- 
ticle size distributions within the regolith. 
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The measurements before and after opposition come 
from different instruments and are too inhomogeneous 
to identify clear trends with phase angle, especially as 
the measurements before opposition (negative phase an- 
gle) have too large error bars. We varied the thermal 
inertia between and 50 J m -2 s -0 5 K _1 to study the 
systematic changes in the data representation of Fig. [5] 
The zero thermal inertia produces a significantly larger 
scatter in the observation-to-model ratios. The T = 50 
calculations showed a very pronounced asymmetry be- 
tween before and after opposition ratios. Our default value 
of 15 J m -2 s~ - 5 K _1 seems to work fine, but somewhat 
smaller or larger values cannot be excluded. The much 
larger scatter in the before opposition ratios is pure co- 
incidence. E.g. the before opposition IRAS data cover a 
much larger range in flux and the most deviating IRAS ra- 
tios on the left side of Fig.|3]belong to the lowest measured 
fluxes. 

The beaming model with p = 0.7 and / = 0.6 influ- 
ences strongly the slope of the observation to model ra- 
tios in the mid-IR (Dotto et al. I2000f) . The flat curves in 
Figs. ^aiid[2Icoiifirm our "default beaming model" . About 
10-20% modifications of the r.m.s. slope value (p) and the 
fraction of surface covered by craters (f) would still be 
in agreement with the PHT-S and the CAM filter mea- 
surements. Miiller (|2002fl demonstrated in different model 
simulations the influences of the thermal parameters with 
wavelength and phase angles. But thermal inertia, beam- 
ing effects and the low level mineralogic emission features 
play together and influence the slope of the thermal emis- 
sion in a complicated way and the different contributions 
are not easy to disentangle. 

4.3. Emissivity analysis 

Barucci et al. 1)1987(1 classified 65 Cybele as CO-type as- 
teroid on basis of visual multi-color and IRAS data. They 
associated this taxonomic type with carbonaceous chon- 
drites. The asteroids 10 Hygiea and 511 Davida, also C0- 
types, were found to show similarities in near-IR and 3 pm 
spectroscopy with carbonaceous chondrite meteorites of 
CI/CM subclass (Hiroi et al. H996|l . Therefore, we concen- 
trated our search for meteoritic analogues to carbonaceous 
chondrites, which represent, chemically and mineralogi- 
cally, the most primitive material in the solar system. 

In Fig. we show the average CVF scan (divided by a 
TPM prediction with e = 1.0) with the error bars derived 
from the two CVF scans {o mean j For comparison, 
we plotted the emissivity of different powdered stony me- 
teorite samples, representing 17 carbonaceous chondrite 
meteorites at small particle sizes of 0-75 pm (Salibury et 
al. I1991al ITMTBI ASTER spectral library 8 ). This library 
contains 60 meteorite samples all at grain sizes of 0-75 pm. 
The sample includes five meteorites of subclass CI, five of 
subclass CM, five of subclass CO and two of subclass C03. 



http:/ /speclib.jpl. nasa.gov 
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Fig. 6. The average CAM-CVF scan was divided by a 
TPM prediction with model emissivity set to 1.0 at all 
wavelengths. The average CVF/TPM ratio is plotted as 
dashed line. The 65 Cybele emissivity curve shows no 
prominent features in the 7 to 15 pm. range. Different car- 
bonaceous chondrites are shown for comparison and ver- 
tically shifted for clarity. 

The emissivity spectra (1 - measured reflectivity) are off- 
set for clarity. 

The 65 Cybele emissivity shows no prominent features, 
but an increase in emissivity between 8.0 and 9.5 pm is 
seen in both down and upward CVF scans and is thus a 
reliable feature, contrary to the part below 7.0 pm where 
the source becomes weak and transient problems hamper 
the measurement (see Sect. I2.3|) . Between 12.5 and 14 pm 
we recognise a very shallow feature which is also present 
in both scans. 

Salibury et al. (1991a) noted that a rather steep in- 
crease in emittance, the so-called Christiansen feature 
(CF), occurs between the wavelength region dominated by 
volume scattering (roughly below 7.5 pm) and the wave- 
length region dominated by surface scattering and rest- 
strahlen bands (roughly beyond 8.5 pm). They placed the 
CF wavelengths of different types of carbonaceous chon- 
drites within the region between about 8.7 and 9.4 pm. 
The CO and C03 types of the carbonaceous (Fig. 
come closest in CF similarity with the 65 Cybele emit- 
tance peak between 9.0 and 9.4 pm, but the data quality 
is not sufficient for a conclusive identification. The two 
C03 stony meteorites in Fig.[l)]are Isna (upper curve) and 
ALHA 77003 (lower curve), the CO meteorites are Lance 
(2x), ALH 83108, Ornans and Warrenton (Salibury et al. 
I1991a[) . The CI and CM subclasses show the CF at slightly 
shorter wavelengths. 

Our results from the emissivity analysis (Sect. I4.2|l 
pointed at predominantly larger particles as cause for the 
far-IR emissivity drop. But the Salibury et al. I|1991afl li- 
brary provides only spectra of small grain sizes. Le Bras & 
Erard 2003 found that the grain size effect can shift the CF 
in wavelengths and change also the contrast. Meteoritic 
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samples with larger grain sizes are therefore needed to 
proove direct connections between the measured spectro- 
scopic features on asteroids and meteoritic samples. 

For the feature between 12.5 and 14 /im we found 
no counterpart yet, but phyllosilicates, like Kaolinite, 
Lepidolite or Serpentine, show emissivity increases at 
around 12.5 /im, similar to what we see here. The fine- 
grained phyllosilicates are one main component of CI/CM 
carbonaceous chondrites, while the CO types are anhy- 
drous materials. This would be an argument in favour of 
similarities with CI or CM types, but again, laboratory 
studies at larger grain sizes are needed to clarify the situ- 
ation. 

Figure El clearly reveals the difficulties of identify- 
ing mineralogic features based on mid-IR spectroscopy. 
Laboratory spectra of different meteorites at various grain 
sizes would be necessary. Concerning the observations of 
asteroids, a high S/N ratio over a large wavelengths range 
is necessary to detect the low level emissivity changes. 
These features would be very difficult to detect from 
ground-based N-band spectroscopy (8 to 13 /im) where at- 
mospheric effects limit the wavelength range and degrade 
the data quality. Additional complications arise from the 
strong slope changes around the thermal continuum en- 
ergy peak at around 14 /im. Poor spectral information to- 
gether with wrongly modelled continuum slopes can then 
easily lead to misinterpretation at the band borders. The 
modelling of the surface temperature distribution, the 
thermal inertia as well as the beaming parameters p and / 
influence the slope of the thermal emission. Our thermal 
values (Sec. 13.111 worked perfectly for the PHT-S and the 
CAM-CVF data sets, but the spectroscopic interpretation 
still remains difficult. 

4.4. Comparison with the STM 

The shape and spin vector information is very important 
for a consistent analysis of the combined observations. 
Although we only used a simple ellipsoidal shape model, 
one can see the shape influence by comparing Fig. with 
Fig. 0] in Sect. 14.11 For Fig. we calculated the model 
values via the STM (Lebofsky et al. I1986[l . using a spher- 
ical shape model. The scatter between data sets taken 
under different aspect angles is much larger. Some mea- 
surements which agree in Fig. 0] are now offset by large 
amounts (e.g. the IRAS 12 and 25 /im observations and the 
CAM filter measurements). Assuming a spherical shape 
(with the given spin vector and thermal behaviour) leads 
to a weighted mean diameter of 280.0 km and an albedo of 
0.047. But now the standard deviations would increase by 
about 60 % for the diameter and about 50 % for the albedo! 
This result clearly demonstrates the importance of using 
spin vector and shape information for combining observa- 
tions from different epochs. Modifying the STM beaming 
value might improve the "Observation/STM" ratios for 
some observations, but different observations would re- 
quire then different STM beaming values. And, for some 
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Fig. 8. The derived diameters and albedos for the CAM 
filter measurements as a function of wavelength, for the 
TPM and also the STM. 



data sets the STM would require a wavelength dependent 
i] to account for the trends in Fig. [7] 

We used the most reliable data set -the CAM filter 
measurements- for a comparison between TPM and STM 
radiometric results. The TPM analysis leads to a con- 
stant diameter and albedo value for all CAM filter mea- 
surements, while the STM produces wavelength depen- 
dent values (Fig. [SJ. The albedo values show the oppo- 
site behaviour. The derived diameters show no trend with 
wavelengths (see also Fig. [21 lower part for a compari- 
son of the observation to TPM ratios over wavelengths) 
and agree within 2-3% (about 6% for the albedo). This 
also confirms our simplified constant emissivity values 
in this wavelength range close to the thermal emission 
peak. Simulations of the 65 Cybele constellation during 
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the CAM filter epoch show that this wavelength depen- 
dent trend continues out to about 30 /im. At longer wave- 
lengths also the STM would produce an almost constant 
diameter. At 30 ^m the STM diameter would be about 
10 % larger. This wavelength dependent effect of the STM 
is mainly produced by the STM correction factor for in- 
frared beaming. Detailed studies of the beaming influence 
with wavelengths and phase angles can be found in Miiller 

rcnna . 

Figures 0] [5] and reveal additional problems of the 
STM. The scatter in the observed fluxes divided by model 
calculations is much larger for the STM (Fig. 01 . The main 
reason for the larger scatter is the assumption of a spheri- 
cal shape. But also the simplified beaming, emissivity and 
phase angle corrections can be recognized in Figs.rj]and0 
The residual beaming effects cause wavelengths dependent 
ratios in Fig.0 easily visible for the IRAS data (★-symbol) 
or the CAM filter data (A-symbol). The offset between 
these 2 sets is caused by the wrong cross section due to the 
spherical shape model of STM. The STM calculation of 
the ratios of Fig.[3]shows a strong asymmetry between the 
ratios before and after opposition. Additionally, the scat- 
ter between the different data sets would be much larger. 
Unfortunately, the overlaying shape effects of the inho- 
mogeneous data did not allow to compare the 0.01mag/° 
phase angle correction factor of the STM with the more 
sophisticated calculations of the true illumination and ob- 
serving geometry by the TPM. 

5. Conclusions 

The recently developed TPM allows to combine ther- 
mal infrared observations with information from visual 
photometry, lightcurves, close-up observations or direct 
measurements. Using 65 Cybele as a "text book exam- 
ple", we applied this powerful TPM in several different 
ways. First, we used re-calibrated observations from litera- 
ture and recently obtained thermal infrared measurements 
to derive the size and albedo of 65 Cybele, taking also 
shape and spin vector information into account. Our de- 
rived effective diameter (273.0±11.9km) and albedo val- 
ues (0.050±0.005) deviate significantly from the IRAS 
based resul ts of 237.3±4.2km and 0.071±0.003 (Tedesco 
et al. 12002(1 . We attribute these differences to the limita- 
tions of the STM (Sects. fO PI and K ty . 

The spectroscopic analysis confirmed the proposed 
"default thermal behaviour" (Miiller et al. I1999[) for large 
main-belt asteroids which is dominated by a surface re- 
golith with very low thermal inertia. The related beam- 
ing model and thermal inertia values are the key ele- 
ments to match even the high accuracy ISO measurements 
at the critical mid-IR range. In return, reliable diameter 
and albedo values can be determined from individual fil- 
ter measurements at any thermal infrared wavelength. No 
band-to-band or wavelength dependent correction factors 
are necessary anymore. Additionally, the thermal emission 
modelling provides the basis for the analysis of mineralogic 



features and allows to investigate relations to meteoritic 
samples. 
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